Absence of kinetic effects in reaction-diffusion processes in scale-free networks 
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We show that the chemical reactions of the model systems of A+A^O and A+B^O when per- 
formed on scale-free networks exhibit drastically different behavior as compared to the same reac- 
tions in normal spaces. The exponents characterizing the density evolution as a function of time are 
considerably higher than 1, implying that both reactions occur at a much faster rate. This is due to 
the fact that the discerning effects of the generation of a depletion zone (A+A) and the segregation 
of the reactants (A+B) do not occur at all as in normal spaces. Instead we observe the formation 
of clusters of A (A+A reaction) and of mixed A and B (A+B reaction) around the hubs of the 
network. Only at the limit of very sparse networks is the usual behavior recovered. 
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The model bimolecular chemical reactions of the A+A 
and A+B type have been heavily studied since the pio- 
neer work of Ovchinnikov and Zeldovich , postulating a 
behavior drastically different than the mean field predic- 
tions, especially in low dimensions and in fractal struc- 
tures. These were verified additionaly by numerical sim- 
ulations, which showed pictorially the dominating effect 
in these two systems, i.e. the generation of the depletion 
zone for the A+A type shown by Torney and McConnell 
, and the spatial segregation of the two types of reac- 
tants in the A+B type shown by Toussaint and Wilczek 
0, In the past twenty years a large number of works 
extended these initial ideas and explained in detail how 
these effects come about, showed the existence of several 
temporal regimes, explained the crossovers between the 
early time and long time behavior, and a wealth of other 
information, rendering these systems as some of the most 
heavily studied systems of interacting particles 0. 

Classically, the mean-field prediction for the density p 
of the surviving particles in A+A and cquimolar A+B 
type reactions is 



P Po 



(1) 



where / = 1, fc is the rate constant, and po is the particle 
density at t — 0. In non-classical kinetics, though, an 
'anomalous' behavior has been observed. In a space with 
dimensionality d the exponent / = d/dc for d < dc, and 
/ = 1 for d > dc- The upper critical dimension dc equals 
2 for A+A, and dc — 4: for A+B. Similarly, anomalous 
behavior has been observed when the diffusion-reaction 
process takes place on different geometries, such as on 
fractals where f = dg/dc, or on dendrimer structures 
1^. In all these cases the limiting value of the exponent 
is / = 1, valid even in infinite dimensions. 

In this Letter we study reactions of the A+A— s-O and 
A+B^O type taking place on scale-free networks. Net- 
works of this type have been recently found to character- 
ize a wide range of systems in nature and society 0, 13 , 
including the Internet, the WWW, chemicals linked via 



chemical reactions, sexual contacts, ecological systems, 
etc. The term scale-free refers to the absence of a charac- 
teristic scale in the connectivity of the nodes comprising 
the network. Thus, each node of the system has k links 
to other nodes in the system with a probability 



P{k) - k- 



(2) 



where 7 is a parameter characteristic of the structure 
of the network. Although these systems are very large 
their diameter is usually small, a property which is usu- 
ally referred to as the 'small- world effect'. The topology 
of such a network is quite complex and leads to a drasti- 
cally different behavior for the above mentioned chemical 
reactions. 

Due to the different structure, as compared to lattices, 
it is of interest to examine a reaction-diffusion process 
on a scale-free network and observe the evolution of the 
two major effects, i.e. generation of the depletion zone 
and segregation of the reactants. Actual situations could 
involve a virus- antivirus reaction of the Internet, or epi- 
demics in a social network. Thus, we performed Monte- 
Carlo calculations of these model systems on a scale-free 
network and monitored the usual parameters, i.e. den- 
sity as a function of time, and the exponent / (equation 
l). 

A scale-free network is constructed as follows. For a 
given 7 value we fix the number of nodes N — 10^ and 
we assign the degree k (number of links) for each node 
by drawing a random number from a power law distribu- 
tion P(k) ^ k^'^ . We do not impose any cutoff value for 
the maximum nodes connectivity, so that k can assume 
any value in the range [1,N]. We then randomly select 
and connect pairs of links between nodes that have not 
yet reached their preassigned connectivity and have not 
already been directly linked to each other. Finally, we 
isolate the largest cluster formed in the network. This 
is identified via a burning algorithm, i.e. we start from 
a randomly chosen node, assign an index i to it, and 
'burn' all of its neighbors, by assigning the same index i 
to them. The process repeats iteratively and when there 
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FIG. 1 : Plots of the reaction progress 1/p— 1/po , as a function 
of time for (a) the reaction and (b) the A+B^O 

reaction on scale-free networks of (left to right) 7 =2.0, 2.5, 
3.0, and 3.5. The initial density was po = 0.5 and po = 0.25 
respectively. All results correspond to networks of N = 10^ 
nodes, except for the dotted lines in (b), where we present 
results for systems oi N = 10* and 10 nodes. The symbols 
represent the simulation data, while continuous lines are the 
asymptotic best-fit lines. 



are no more neighbors to be burnt in this cluster, a new 
non-burnt node is chosen again, an index i+1 is assigned 
to it, etc, until all nodes have been visited with this pro- 
cess. The largest cluster can now be easily identified by 
counting the number of nodes with the same index. 

The reaction mechanism is similar to the one routinely 

used in the literature for the case of regular lattices. For 
the A-hA^O reaction, an initial fraction of the nodes po 
of the largest cluster is randomly chosen. These nodes 
arc occupied by A particles. One particle is chosen at 
random and the direction of its move is also chosen at 
random with equal probability to one of its neighbor sites 
(nodes directly linked to its present position). The total 
time advances by 1/A^. If this node is already occupied 
by another particle the two A particles are annihilated, 
otherwise the particle is moved to the new node. For the 
A-|-B^O reaction, an equal initial fraction po of A and B 
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FIG. 2: Asymptotic reaction rates (slopes of the lines in 
figure 1), as a function of 7, for the A+A^O reaction (open 
symbols) and the A^-B— >0 reaction (filled symbols). 



particles are randomly placed on the network. Particles 
move using the same algorithm described above, but now 
when a particle tries to move to a node where a particle of 
the same type resides the move is not allowed, although 
the time advances. When an A particle encounters a B 
particle these two particles annihilate. We monitor the 
population of particles as a function of time. 

In figure 1 we present the evolution of the particle con- 
centration as a function of time for the (a) A-l-A and (b) 
A-l-B reactions. Results for networks with different 7 
values are presented, where 7 ranges from 2 to 3.5, rep- 
resenting a varying degree of nodes connectivities (from 
dense to sparse networks). As we can see, the curves 
follow power-law behavior with two distinct regimes. In 
all cases, there exists a crossover between the early-time 
regime and the asymptotic limit. The location of the 
crossover point increases with 7. For 7 = 2 the crossover 
takes place as early as 10 steps or less, while for 7 = 3.5 
it is of the order of 1000 steps. The reaction rate, also, 
is much faster for lower 7 values, with the concentration 
falling to 10~^ in only 50 steps for 7 = 2, for both types 
of reaction. In figure lb and for 7 = 3 we also give results 
for smaller networks (A'" = 10^ and A'^ = 10^), in order 
to ascertain that there are no finite-size eff'ects for the 
calculation of the asymptotic slopes. 

These observations are in contrast with reaction 
schemes on regular lattices. This can be clearly seen 
in the calculation of the asymptotic slopes for the curves 
of figure 1. These slopes (corresponding to the exponent 
/ of equation 1) are presented in figure 2, as a function 
of 7. There is a small difference between the slopes of 
the A+A and the A+B types (with the latter slightly 
larger), but in all cases the slope is greater than 1. The 
exponent / can acquire very large values (/ = 2.75 for 
7 = 2) and is a monotonically decreasing function of 7. 
As already stated, on all other geometries studied in the 
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FIG. 3; (a) A+A reaction: Percentage Qaa of contacts be- 
tween A particles over the total number of possible contacts 
Naa/N{N — 1) as a function of time, (b) A+B reaction: Per- 
centage of AB contacts over (AA~f BB) contacts as a function 
of time. The line at Qab ~ 1 corresponds to complete mixing. 
The 7 values are as marked. 



literature, the exponent / is lower than 1, while on scale- 
free networks we see that the value / = 1 is reached only 
asymptotically for networks with large 7. 

In order to understand the spatial distribution of the 
particles, we calculate the number of close contacts they 
form (we consider a close contact as the existence of a 
link between two network nodes). In the case of the 
A+A reaction, for fixed time, we measure the fraction 
Qaa defined as the number of contacts iV^^ between A 
particles over the total possible number of contacts, i.e. 



Qaa — 



Naa 



N{N ~ 1) 



(3) 



The value of Qaa = 1 corresponds to the extreme case of 
all particles forming one cluster, while a decrease of this 
value indicates the existence of a depletion zone (parti- 
cles are placed apart from each other). In figure 3a we 
calculate the fraction Qaa for different networks as a 
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FIG. 4: Lattice analogue pictorial of the particles configura- 
tion in space for the A+B-^0 reaction. Notice that particles 
of different kind gather in the same region, in the close vicinity 
of the hubs. 



fimction of time. When 7 = 2 or 2.5 a depletion zone 
initially forms, but immediately the number of contacts 
increases with time. Most of the particles are clustered 
in a small region and as time advances this clustering 
increases, i.e. particles continue to gather in positions 
close to each other, in the vicinity of the most connected 
nodes (hubs). This is in striking contrast to the well- 
established formation of a depletion zone on regular lat- 
tices. For 7 = 3 we can see that the number of contacts 
remains almost constant with time, while for 7 = 3.5 we 
can recover the formation of the usual depletion zone, 
where the number of contacts decreases monotonically 
with time. 

Similar conclusions are found for the A+B reaction 
(figure 3b). In this case we measure the number of con- 
tacts between unlike particles compared to the number 
of contacts between particles of the same type. We use 
the ratio Qab defined as 



Nab 



{AB 



Naa + N, 



BB 



(4) 



The value of Qab tends to when segregation of like 
species is formed, while it tends to a value of 1 when 
complete mixing of A and B occurs. When 7 = 2 there 
is clearly no segregation of like particles, even at early 
times, but both A and B particles mix together. For 
7 = 2.5 and 7 = 3.0 there are hints of segregation, since 
the number of unlike particles contacts decreases com- 
pared to like particle contacts, but soon particles start to 
aggregate. The time needed for the transition to mixing 
and the extend of the depletion zone increase with 7. At 
7 — 3.5 the classical picture with a clear segregation of 
particles reappears, where like particles gather in clusters 
with a scarce presence of unlike particles at all times. 

This peculiar behavior can be attributed to the struc- 
tural characteristics of the scale-free networks. Such net- 
works are known to heavily rely on the existence of hubs. 
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which are nodes with a large number of connections. 
Moreover, such networks have been shown to have an ex- 
tremely small diameter, of the order In(lniV) These 
two factors can explain the results of the present work. 
The small diameter of the network causes the majority of 
the particles to be at a close distance to each other, and 
most of the nodes can reach a hub through a small num- 
ber of links. Thus, we expect that diffusion will very soon 
bring the particles close to each other and they will react 
within a short period of time. This is indeed the case of 
figure 1, where we can see that the process ends rapidly. 
The importance of the hubs diminishes as we increase 7, 
since now nodes are less connected to each other, and 
the familiar notions of depletion zones and segregation 
are restored for 7 — 3.5. One can picture this process as 
a biased walk of the particles towards the hubs, where 
the probability of encounter of another particle is much 
higher than in the case of a regular lattice walk. 

These concepts correspond to the idea of continuous 
mixing. For the A-l-A reaction and low 7 values there 
exists no depletion zone. This is due to the fact that 
neighboring A particles annihilate initially in a very fast 
rate, but as time advances the diffusion of the A particles 
brings them close to the hubs, which dominate the struc- 
ture. Thus, there are always pairs of A close to each other 
(absence of depletion). For larger 7 values, i.e. when the 
networks are more sparse, the connectivity of the hubs 
diminishes, and the regular result of a depletion zone is 
reproduced. 

In figure 4 we present a lattice analogue pictorial of 
the A-l-B case. Although unlike particles that were close 
at t = have been annihilated, at longer times there is 
still no segregation, since all particles, independently of 
their type, move towards the hubs. We thus have the 
formation of clusters made of A and B particles, which 
are located close to the hubs. This means that A and 
B particles are always mixed because of the underlying 
structure and the reaction advances at a constant rapid 
rate. We have indeed verified in our simulations that the 
majority of the annihilation reactions take place on the 
hubs. Notice also that diffusion on scale-free networks 



is non-Markovian, since the walk depends on the exact 
particle location, due to the largely varying connectivity 
of the nodes. 

A similar situation is encountered in another system, 
that of Levy walks, where it was observed that the rare 
long-distance jumps break the formation of the segre- 
gated regions, acting in effect like a stirring mechanism 

Summarizing, we have presented evidence that 
reaction-diffusion processes in scale- free networks are dif- 
ferent in their nature compared to lattice models, ex- 
hibiting rapid reaction rates (/ > 1). This is due to the 
small diameter of the networks and the existence of the 
hubs. These differences are more pronounced in compact 
networks of low 7 values, while for sparse networks, e.g. 
7 > 3.5, the behavior is the same as for regular lattices. 
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